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FOREWORD 


Technical  Report  No.  2303-1  was  prepared  by  the  Stanford 
Electronics  Laboratories,  Stanford  University,  on  Air  Force  Contract 
AF04(695)-305  under  Task  No.  318201  of  Project  No.  3182,  "University 
Program  for  Vehicle  Detection  and  Defense."  The  work  was  administered 
under  the  direction  of  the  Air  Force  Space  Systems  Division,  Air  Force 
Systems  Command.  Capt .  Robert  Eaglet  was  Project  Engineer  for  the 
Division. 

The  studies  presented  covered  the  period  from  July  through  November 
1963.  The  research  activity  was  conducted  by  Joseph  W.  Goodman  at  the 
Systems  Techniques  Laboratory  at  Stanford. 

This  is  an  interim  technical  report  on  one  phase  of  the  work  under 
Contract  AF04(695)-305. 
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ABSTRACT 


When  laser  light  strikes  a  diffuse  object,  such  as  paper,  the 
scattered  light  has  been  observed  to  possess  a  granular  spatial  structure. 
The  statistical  properties  of  these  so-called  "sparkle  patterns,"  as  seen 
by  an  observer  in  the  far  field  of  the  scattering  spot,  are  investigated. 

The  first  order  statistics  of  the  observed  electric-field  strength, 
the  observed  light  intensity,  and  the  observed  light  phase  are  examined. 
The  electric  field  is  reasoned  to  be  a  complex  normal  random  variable; 
the  intensity  a  real,  exponentially  distributed  random  variable;  and  the 
phase  a  uniformly  distributed  random  variable.  Higher  order  statistics 
of  these  random  processes  are  also  discussed.  The  autocorrelation  func¬ 
tions  of  the  complex  field  and  the  intensity  processes  are  investigated, 
and  that  of  the  electric  field  is  found  to  be  proportional  to  the  Fourier 
transform  of  the  light-intensity  distribution  incident  on  the  scattering 
surface. 

Spatial  averages  of  the  light  intensity  are  considered  and  are  found 
to  converge  to  corresponding  ensemble  averages  when  either  the  area  of 
the  scattering  spot  or  the  averaging  area  grows  large. 

The  influence  of  light  spectral  content  on  sparkle  patterns  is 
investigated  in  some  detail.  The  cross  correlation  of  patterns  produced 
by  different  monochromatic  light  components  is  evaluated.  The  degree  of 
correlation  is  found  to  depend  on  the  frequency  difference,  the  roughness 
of  the  scattering  surface,  the  angle  of  incidence  of  the  laser  beam,  and 
the  position  of  the  observer.  Light  spectral  components  of  nonzero  width 
are  considered,  and  the  influence  of  bandwidth  on  the  ability  of  an 
observer  to  distinguish  a  sparkle  pattern  is  examined,  it  is  found  that 
sparkle  patterns  exist,  for  light  of  any  bandwidth,  but  t f mo-invariant 
patterns  can  be  seen  only  when  the  bandwidth  is  less  than  a  -ertain  limit 
which  depends  on  the  roughness  of  the  scattering  surface,  the  angle  of 
incidence  of  the  laser  beam,  and  the  position  of  the  observer. 

Publication  cf  this  technical  documentary  report  does  not  constitute 
Air  Force  approval  of  the  report's  findings  or  conclusions.  It  is 
published  only  for  the  exchange  and  stimulation  of  ideas. 
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sparkle  on  laser  ranging  ami  vel set ty-measuroment  systems  will  be  inves 
tigated. 
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th 

where  DQ  is  the  distance  from  the  k  scatterer  to  the  observation 
point, 

2  th 

| E. |  is  the  intensity  of  the  light  incident  on  the  k  scatterer, 

—  and 

i  1 2 

|ek|  is  the  intensity  of  that  portion  of  the  light  at  the  observa- 

—  tion  point  which  arises  from  the  k**1  scatterer. 

The  scattering  cross  section  is,  in  general,  a  function  of  both  the 
light  frequency  and  the  observation  point.  It  is  assumed  that  the 
observation  region  is  sufficiently  small  and  the  distance  DQ  suffi¬ 
ciently  large  so  that  is  constant  within  that  region  for  any  given 

frequency.  By  using  the  concept  of  cross  section,  the  scattering  ele¬ 
ments  are  considered  to  be  isotropic  point  radiators.  Further  assumptions 
on  the  statistical  properties  of  the  elementary  scatterers  will  be  dis¬ 
cussed  below. 

At  the  scattering  surface  an  "effective"  scattered  field  distribu¬ 
tion  is  defined  by  projecting  the  point  radiators  onto  a  plane  X-Y  at 
the  surf  ace,  which  is  parallel  to  the  observation  plane  U-V.  A  Z  axis 
is  defined  running  perpendicularly  between  the  origin  of  the  X-Y  plane 
anc  the  origin  of  the  U-V  plane,  as  shown  in  Fig.  1.  The  point  z  =  0 
is  de'ined  to  coincide  with  the  origin  in  the  X-Y  plane,  while  the 
origin  of  the  U-V  plane  is  assumed  to  lie  at  z  =  D,  w  e  D  repre¬ 
sents  the  distance  between  the  two  planes. 

An  effective  complex  field  distribution  for  the  projected  scatterers 
in  the  X-Y  p’ ane  may  be  written  as 

K 

E(x,  y;  t)  exp  (j&k)  6(x  -  xR,  y  -  yfc)  exp  (jasQt)  (2.2) 

k=l 
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FIG.  1.  THE  GEOMETRY  OF  THE  SCATTERING  AND  OBSERVATION 
PLANES. 


where  K  is  the  number  of  scatterers  contributing  to  the  observed 
field  in  the  U-V  plane; 


ak  is  equal  to  ysk/4i;  jEj  and  is  a  measure  of  the  strength 
of  the  kth  scatterer; 


th 

Pk  is  the  phase  of  the  radiation  from  the  k 
measured  at  its  projected  point  in  the  X-Y 


scatterer,  as 
plane; 


(xk’yk) 


th 

is  the  location  of  the  projection  of  the  k  scatterer  in 
the  X-Y  plane; 


u>0  is  the  angular  frequency  of  the  electromagnetic  radiation; 
and 


5(  •  ,  *  )  is  a  two-dimensional  Dirac  delta  function. 
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To  distinguish  between  the  spatial  and  time  variations  of  the  field, 
Eq.  (2.2)  can  be  written 


E(x,  y;  t)  =  E0  (x,  y)  exp  (ju>0t)  (2.3) 


from  which  it  is  evident  that 


E0  y ^  lEi  ^xk'  yk^l  exp  '  Xk'  y  "  yk^  (2‘4) 

—  k=l  — 


In  the  treatment  which  follows,  the  function  EQ  (x,  y)  is  con¬ 
sidered  to  be  a  random  process  over  an  ensemble  of  possible  scattering 
distributions  in  the  X-Y  plane.  The  following  assumptions  are  made: 

1.  The  scatterers  are  randomly  distributed  over  the  X-Y  plane  with 
uniform  probability  (as  mentioned  earlier). 

2.  The  scattering  cross  section  of  a  given  scatterer  is  a  random 
variable,  and  all  scatterers  have  statistically  independent  and 
identically  distributed  cross  sections. 

3.  Since  nearly  all  surfaces  are  extremely  "rough"  compared  to  a 
light  wavelength,*  the  phase  angles  6^  are  uniformly  distributed 

in  probability  over  the  interval  (0,  211 ) . 

4.  Each  Pk  is  statistically  independent  of  all  others. 

The  final  assumption  of  statistical  independence  might  be  improved 
by  allowing  correlation  to  exist  between  the  £  of  neighboring  scat- 
terers.  However,  providing  that  the  correlation  exists  over  an  interval 
which  is  small  compared  to  the  total  extent  of  the  scattering  "pot,  no 
major  changes  in  the  results  presented  here  are  anticipated. 

'*tf  the  scattering  surface  is  not  rough  compared  to  a  wavelength,  the 
results  presented  here  apply  only  to  that  portion  (possibly  small)  of 
the  observed  radiation  which  is  not  specular  reflection.  However,  the 
results  do  apply  to  quasi-specular  reflection  from  metallic  objects 
caused  by  reflection  from  many  small,  optically  smooth  areas,  each  cf 
which  is  randomly  distributed  in  depth  over  several  wavelengths. 
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III.  FIRST-ORDER  STATISTICS  OF  THE  SCATTERED  LIGHT 
IN  THE  OBSERVATION  PLANE 

A.  AN  EXPRESSION  FOR  THE  COMPLEX  FIELD  STRENGTH  IN  THE  OBSERVATION 
PLANE 

The  relationship  between  the  electric  field  in  the  X-Y  plane  at 
the  scattering  surface  and  the  electric  field  in  the  observation  region 
of  the  U-V  plane  is  investigated  in  Appendix  A.  It  is  shown  there 
that  a  Fourier  transform  relationship  exists  between  the  electric  fields 
in  the  two  planes;  specifically,  if  s(u,  v;  t)  represents  the  complex 
electric  field  strength  in  the  U-V  plane,  then 


I  I  Eo(x'  exp  [j-1  ^UX  +  Vy^jdx  dy 


provided  the  following  conditions  are  met: 


1.  x/D  «  1  and  y/D  «  1  for  (x,  y)  in  the  scattering 

spot. 

2.  u/D  «  1  and  v/D  «  I  for  (u,  v)  in  the  observation 

region. 

2  2 

3.  2Hx  /AD  «  1  and  2Hy  /AD  «  1  (i.e.,  the  observation  region  is 


in  the  far  field  of  the  scattering 
spot). 


4.*  flu^AlAD3  «  1  and  Hv4/4AD3  «  1. 
For  simplicity  Eq.  (3,l)  is  rewritten: 


r  /  2  2n 

e(u,  v;  t)  =  a  exp  (ju>0t)  exp  |^-j  -  ^je0(u.  v 


)  (3.2) 


Note  that  this  assumption  does  not  require  the  scattering  spot  to  be  in 
the  far  field  of  the  receiving  optics  (observation  region). 
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where 


e0(u,  v)  =  J  J  E0(x,  y)  exp  jjj  ||  (ux  +  vy)j  dx  dy  (3.3) 


and 


exp 


(  .  21®\ 

V3  T) 


jAD 


(3.4) 


The  random  parameters  of  the  scattering  surface  become  explicitly 
evident  when  (2.4)  is  substituted  into  (3.3),  with  the  result 

K  t _ 

e0<u'  v)  !fi(V  yk}l  exp 

k=l 


Pk  + 


2):  , 

to  K  + 


vyk) 


} 


(3.5) 


From  this  expression  and  the  assumptions  on  the  statistical  proper¬ 
ties  of  the  scattering  surface,  it  is  now  possible  to  derive  the  relevant 
statistical  properties  of  the  scattered  light. 

B.  STATISTICS  OF  THE  COMPLEX  FIELD  STRENGTH 

From  the  expression  (3.5)  it  is  evident  that  the  complex  electric 
field  eQ(u,  v)  is  given  by  a  sum  of  randomly  phased,  random-amplitude 
phasors.  Because  of  the  statistical  independence  and  uniform  distribu¬ 
tion  of  the  random  phases  the  following  assertion  can  be  proved 

(Ref.  5,  p.  362]: 

As  the  number  K  of  elementary  scatterers  in  the  scattering  spot 
increases  without  bound,  the  random  process  s^(u,  v)  becomes  a 

complex  normal  process.  The  real  and  imaginary  parts  of  this  process 
are  identically  distributed  with  zero  means  and  identical  variances 

2 

which  will  be  designated  cr  /2.  At  any  single  point  (u,  v),  the 
real  and  imaginary  parts  of  the  process  are  statistically  independent. 
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It  is  assumed  throughout  that  the  number  of  scattering  elements 
within  the  spot  is  indeed  so  large  that  €Q(u,  v)  may  be  considered  to 

have  normal  statistics.  Thus,  if  fQ(u,  v)  is  written  as  the  sum  of 
real  and  imaginary  parts, 

e0(u,  v)  =  e^u,  v)  +  j  e2(u,  v)  (3.6) 

the  „yint  probability  density  function  of  € ^  and  way  be  written 

as 


e2> 


1 

— ~  exp 
h<t 


(3.7) 


From  the  relationship  (3.2)  between  e(u,  v;  t)  and  €0(u,  v),  it 
can  be  seen  that  for  a  fixed  (u,  v;  t)  the  former  is  also  a  normal 
process.  However,  as  will  be  seen  in  Chapter  IV,  at  any  fixed  instant 
of  time,  e(u,  v;  t)  in  general  will  have  statistics  which  are  non¬ 
stationary  over  the  U-V  plane,  while  £q(u.  v)  is  a  stationary 
process. 


C.  STATISTICS  OF  THE  LIGHT  INTENSITY  AND  PHASE 

The  light  intensity  n  the  observation  plane  is  defined  as  the 
square  of  the  modulus  of  '-he  complex  electric  field  strength, 

I(u,  v)  £  |e(u,  v;  t)|2  (3.8) 


From  (3.2)  it  is  clear  that 


The  intensity  is  an  important  physical  quantity  to  consider,  in  that 
the  majority  of  light  detectors  have  a  response  that  is  linearly  pro¬ 
portional  to  intensity. 
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Although  most  detectors  are  Insensitive  to  the  phase  of  the  received 
radiation,  applications  do  exist  in  which  the  phase  plays  an  important 
role  (e.g.,  heterodyne  receivers).  Accordingly,  the  phase  of  the 
received  radiation  is  defined  by 


v) 


-1 

tan 


e2(’-i,  v) 
^(u,  v) 


(3.10) 


and  its  statistical  properties,  as  well  as  those  of  the  intensity,  will 
be  investigated. 

To  find  the  first-order  statistics  of  I  and  y,  the  transformation 
defined  by  (3.9)  and  (3.10)  must  be  applied  to  the  random  variables 
and  Sg.  The  inverse  of  this  transformation  is  easily  seen  to  be 


r1/2 


6.  =  -  cos  -Jr 

laj 


I1/2 


=  -  sin  'If 

‘  1*1 


(3.11) 


from  which  the  Jacobian  is  found  to  be 


|J|  = 


d  I 

d  >Jr 


5  I 


Se2 

5T 


i 

2|a| 


(3.12) 


The  joint  probability  density  function  of  I  and  if  may  now  be  found 
from  the  relation  [Ref.  6,  p.  38] 


■...‘■•♦’-H  V,(g 


I1/2  I1/2  \  /  . 

cos  ijr,  -  sin  -Jr  (3.13) 


111 
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The  result  is 


tf)  = 


1 

2H[a|V 


exp 


(3.14) 


Integrating  over  the  appropriate  variables  reveals  the  marginal  densities 
of  the  intensity  and  phase  to  be 


and 


f 

1 

(  I  \ 

ISO 

1  I2  2  exp 

1 2.1  cr 

( 

\  |a|V/ 

0 

otherwise 

c 

i 

2n 

0  g  *  S  211 

0 

otherwise 

V 

These  results  may  be  summarized  as  follows: 


(3.15) 


(3.16) 


The  light  intensity  l(u,  v)  is  a  real,  exponentially  distributed 
random  variable,  while  the  light  phase  \Jr(u,  v)  is  a  uniformly 
distributed  random  variable  on  the  interval  (o,  2!l). 


Note  that 


pi,*(l-  =  pi(l)  P*W 


(3.17) 


from  which  it  follows  that: 

the  intensity  and  phase  at  the  point  (u,  v)  are  statistically 
independent  random  variables. 

The  moments  of  the  intensity  can  be  shown,  from  (3.15),  to  be 

!V  =  |a|2V(r2v  r(v  +  l)  (3.18) 
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;V.  HIGHER  ORDER  STATISTICS  OF  THE  SCATTERED  LIGHT 


A.  INTENSITY-DISTRIBUTION  FUNCTIONS  AND  CORRELATION  FUNCTIONS 

In  the  previous  chapter  the  firsl-order  statistics  of  the  complex 
field  strength,  the  intensity,  and  the  phase  of  the  scattered  light  were 
found.  A  complete  statistical  description  of  these  quantities  at  any 
single  point  in  space  is  known  therefore,  but  not  the  statistical  rela- 
-xonships  between  their  values  at  two,  three,  or  N  different  points 
in  space  (i.e.,  the  "higher  order"  statistics).  It  is  a  remarkable 
property  of  the  normal  random  process  that  once  the  so-called  "auto¬ 
correlation  function"  of  the  process  is  known,  the  statistics  of  all 
orders  are  completely  determined.  In  order  to  more  fully  describe  the 
statistics  of  sparkle  patterns,  the  concept  of  autocorrelation  function 
and  the  related  concept  of  "intensity-distribution  function"  are  intro¬ 
duced. 

The  two-dimensional  intensity-distribution  function,  S(x,  y),  of  the 
scattered  light  intensity  at  the  point  (x,  y)  is  defined  as  the 
ensemble  average  of  the  scattered  light  intensity  at  that  point.  It  is 
assumed  that  the  elementary  scatterers  are  so  dens,  ly  packed  that  S(x,  y) 
may  be  considered  to  be  a  continuous  function,  in  which  case  it  can  be 
seen  that 


s(x,  y)  =  jjf  |e1(xi  y)|2  (4.1) 

where  s  is  the  ensemble  average  of  the  scattering  cross  section  of  the 
elementary  radiators.  Thus  the  functional  form  of  S(x,  y)  is  identical, 
up  to  a  multiplicative  constant,  with  the  functional  form  of  the  intensity 
distribution  across  the  scattering  spot.  The  intensity-dictribution  func¬ 
tion  has  the  particular  property  that,  when  integrated  over  the  entire 

2 

X-Y  plane,  it  yields  the  constant  cr  which  is  twice  the  variance  of 
the  normal  distributions  discussed  earlier. 
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The  two-dimensional  autocorrelation  function  of  an  arbitrary  complex¬ 
valued  random  process  £(u,  v)  is  defined  by 

F^(Su,  8v)  =  £(u,  v)  i*(u  -  bu,  v  -  5v)  (4.2) 

where  the  superbar  indicates  an  ensemble  average.  A  process  is  said  to 
be  wide  sense  stationary  if  the  correlation  function  depends  only  on 
(5u,  8v);  while  if  it  depends  on  (u,  v)  as  well  as  (8u,  5v),  the 
process  is  nonstationary. 

The  following  important  properties  of  any  autocorrelation  function 
should  be  noted: 


1.  (6u,  8v)  is  in  general  complex. 

2.  | ft.  (8u,  8v)|  C  |ft.  (0,  0)|  . 

— i 

3.  ^  °)  is  real  and  equals  twice  the  average  intensity  of  the 

process  £( u,  v). 


Using  the  definition  (4.2),  the  autocorrelation  functions  of  the 
random  processes  eQ(u,  v),  e(u,  v;  t),  and  l(u,  v)  will  now  be 
investigated.  From  (3. 5)  and  (4.2)  the  autocorrelation  function  of 
e0(u,  v)  is 


P'g  <6u’  8v)  =2^  lfi(V  exp 

—  k=l 

Note  that  this  autocorrelation  function  does  not  depend  on  (u,  v),  so 
the  random  process  eQ(u,  v)  is  wide-ser.se  stationary.  Further  note 
that  if  the  scatterers  are  sufficiently  dense  in  the  scattering  spot, 
the  sum  may  be  approximated  by  an  integral,  with  the  result 

(5u,  6v)  =  J  J  ^  |  EjCx,  y)|  2  exp  j^j  ||  (x&u  +  y&v)j  dx  dy  (4.4) 


j 

— 


211  / 

DA  W 


,  &u 


+  uk5v) 


(4.3) 
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Recalling  the  definition  (4.l)  of  the  intensity-distribution  function, 
it  follows  that: 

The  intensity-distribution  function  and  the  autocorrelation  func¬ 
tion  of  the  random  process  e0(u,  v)  are  a  Fourier  transform  pair; 
specifically  — 


(Su,  &v)  =  J  J"  S(x,  y)  exp  j  gjj  (x6u  +  y5v)J  dx  dy  (4.5) 


Since  the  intensity-distribution  function  is  proportional  to  the  intensity 
of  the  incident  light  at  the  scattering  surface,  the  functional  form  of 
the  autocorrelation  function  can  be  found  directly  once  the  intensity 
distribution  across  the  scattering  spot  is  known.  Note  that  for  a  given 
intensity-distribution  function,  the  width  of  the  corresponding  auto¬ 
correlation  function  depends  on  the  wavelength  A.  The  longer  the  wave¬ 
length,  the  wider  the  peak  of  the  correlation  function. 

The  autocorrelation  function  of  e(u,  v;  t)  is  now  examined 
briefly.  Equations  (3.2)  and  (4.2)  yield 

P..(6u,  Sv)  =  |  a| 4  F.^  (6u,  5v)  exp  j"-  j  (2u&u  +  2v6v  -  5u2  -  5v2)J 

(4.6) 

Since  this  expression  is  a  function  of  (u,  v)  as  well  as  (&u,  5v), 
the  process  e(u,  v;  t)  is  nonstationary.  However,  it  can  be  shown 
that  if  the  scattering  spot  is  in  the  far  field  of  the  observation 
region,  then  the  dependence  on  (u,  v)  is  eliminated  and  the  process 
is  wide-sense  stationary. 

Turning  now  to  the  autocorrelation  function  of  the  intensity  random 
process,  from  (3.9)  and  (4.2)  it  is  seen  that 

Rj(8u.  Sv}  =  I  a|  4  |e0(u,  v)|2  |e0(u  -  5u,  v  -  6v|"  (4.7) 
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Expanding  £g(u,  v)  in  its  real  and  imaginary  parts 


£0(u>  v)  =  ei(u>  v)  +  J  «2(u.  v) 


(4.8) 


yieius 


pj(8u,  6v)  =  ja|4  e^(u,  v)  e^(u  -  8a,  v  -  8v) 


+  e^u,  v)  e^u  -  Su,  v  -  8v) 


2  2 

+  e2(u,  v)  €i^u  '  ^u’  v  " 


+  <^(u,  v)  £g( u  -  8u,  v  -  8v) 


(4.9) 


Now  if  x^,  xQ,  x3,  and  x^  are  real,  zero-mean,  normal  random  vari¬ 
ables,  then  the  following  relation  is  known  to  be  true  [Ref.  5,  p.  343): 


X1X2X3X4  =  X1X2  X3X4  *  X1X3  X2X4  +  X1X4  X2X3  (4-10) 


Using  this  fact,  it  is  not  difficult  to  reduce  (4.9)  to 

ft 


^Su,  5v)  =  |  a |  4  (T4  +  |ft_  (8u,  Sv)|2j  (4.1l) 


which  provides  a  simple  relation  between  the  autocorrelation  functions 
of  eQ(u,  v)  and  l(u,  v).  Note  that  the  random  process  l(u,  v)  is 
wide-sense  stationary. 

Three  important  intensity-distribution  functions  and  their  corre¬ 
sponding  autocorrelation  functions  are  listed  below: 
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3.  Circular  Spot 


s(x,  y) 


P  (8u,  8v) 
€0 


rff2  2  2  2 

— r  0  <  x  +  y  £  p 

np2 


otherwise 


ji(^  V^>2  +  (Sy)2) 

^  ^(5u)2  +  (8v)2 


2(7 


(4.14) 


Finally,  some  quantities  are  defined  which  will  be  of  considerable 
use  in  later  chapters. 

The  cross-correlation  function  of  two  complex-valued  random  proc¬ 
esses  ^(u,  v)  and  §2(u.  v)  is  defined  by 


'H  .  (Su,  8v)  =  §  (u,  v)  §*(u  -  8u,  v  -  8v )  (4.15) 

s12  —  — 


In  words,  the  cross-correlation  function  provides  a  measure  of  the 
statistical  similarity  of  the  first  random  process  g1  at  the  point 
(u,  v)  and  the  second  random  process  at  the  point  (u  -  8u,  v  -  8v). 

A  cross-correlation  function  need  not  have  its  maximum  value  at 
8u  =  8v  =  0, 

The  cross-intensi ty-distribution  function  is  defined  as  the  two- 
dimensional  Fourier  transform  of  the  cross-correlation  function, 


sf  (*. 
S12 


y) = r/v6u’ 6v)  ex?  t  ^ (x5u + y8v)]  ^ 


(4.16) 
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The  cross-intensity-distribution  function  is  in  general  a  complex 
quantity. 


B.  HIGHER  ORDER  STATISTICS  OF  THE  COMPLEX  FIELD  STRENGTH 

In  this  section  the  joint  statistics  of  the  real  and  imaginary  parts 
of  the  complex  field  strength  at  A  different  points  in  space  are 

examined.  Since  these  quantities  have  been  seen  to  be  normally  distrib- 
..ted ,  the  statistical  properties  of  the  field  at  N  different  points 
in  space  are  completely  described  by  a  2N-dimensional  normal  distribution, 


pn(c 


11-  12’  21’  22’ 


ht^) 


(2n)N  |a|x/2 


(4.17) 


where 

th 

is  the  real  part  of  the  field  strength  at  the  k  point; 

th 

ek2  is  the  imaginary  part  of  the  field  strength  at  the  k  point; 
€  is  a  column  matrix  of  the  e,  defined  by 


'll 
'  12 
l21 
'22 


N1 


N2 


(4.18) 
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is  the  transpose  of  e; 


A  is  the  covariance  matrix  of  the  e,  defined  by 


A  =  e  = 
- - 1 


sll€ll 

eHei2 

.  .  . 

£lieNl 

'11"N2 

el2C21 

G12ei2 

£  £ 

12W1 

£  £ 

12  K2 

C  fi 

Will 

c  £ 

XI  12 

•  •  • 

'SieNl 

€N1€N2 

S'2ell 

r  € 

X2  i2 

.  .  . 

€N2£N1 

S’2SN2 

(4.19) 


A  ^  is  the  inverse  of  A;  and 
j A j  is  the  determinant  of  A. 

The  elements  of  the  correlation  matrix  can  be  found  from  the  auto¬ 
correlation  function  of  the  process  e0(u,  v)  by  means  of  the  following 
relationship,  which  is  not  difficult  to  prove; 


Rc0(uk  -  V  vk  ■  vi)  =  fo(V  vk}  fo(ur  VP 


"  £kl€£l  +  ek2£|!2^j  +  j  e£iek2  ‘  £^2€kl  (4-20) 


"  2  eklCn  *  J  2  e£l&k2  =  2  6k2Si2  '  j  2  ef!2-kl 


Note  that  the  imaginary  part  of  the  autocorrelation  function  is  propor¬ 
tional  to  the  cross  correlation  between  the  real  and  imaginary  parts  of 
the  field  strength  at  the  two  different  points  involved.  It  was  pre¬ 
viously  shown  that  the  real  and  imaginary  parts  of  the  field  strength 
at  any  one  point  are  uncorrelated,  but  it  is  evident  from  (4.20)  that 
the  real  ar.d  imaginary  parts  at  different  points  may  be  correlated.  If, 
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however,  the  laser  spot  is  symmetrical  about  its  center,  the  autocor¬ 
relation  function  is  entirely  real,  and  the  real  and  imaginary  parts  of 
the  field  strength  are  statistically  independent  processes. 
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The  Jacobian  of  the  transformation  can  be  shown  to  be 


M  =  TTj4 
4ial 


(4.23) 


and,  using  the  higher  order  anaiogy  of  (3.13), 


exp 


PllV'l*2  *•  2'  V 


(4.24) 


where  in  this  case 


€  - 


111 


T1/2 

Ij'  COS 


I V2 


sin  * 


T1/2 

1 2  cos 


l\/2  sin 


(4.25) 


It  is  convenient  at  this  point  to  have  separate  symbols  for  the 
real  and  imaginary  parts  of  R  (5u,  6v).  Accordingly, 


R  (5u,  &v)  =  R  (6u,  6v)  +  j  R„(Su,  6v) 
e0  1  i 


(4.26) 
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The  covariance  matrix  can  now  be  written 


(4.27 


The  determinant  |A|  is 


l&l 


ana  the  inverse  matrix 


can  be  shown  to  be 


(4.28 


2 

a 


0 


-R 


1 


-R 


2 


0 


-R 


1 


2 

cr 


2 

cr 


-R 


1 


0 


(4.29 


R 


2 


R 


1 


0 


2 

cr 


Carrying  through  the  computation  required  by  (4.24)  yields 


VaW1''  '2' 


■<l2lz  l-il 


TJi  (li 


<■» 

TJi 


2|“I2  l£l 


‘I*  *2* 


>6'2  i*i4  Id 1/2 


*2  ’  *i>  ♦  r2  sln  ('  2  ‘  ■' 2 ) ), 


(4.30 
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Appropriate  integrations  over  (•■>•  ,  *  )  and  (i  ,  I„)  yield  the  desired 

X  a  x  ^ 

marginal  densities  [cf.  Ref.  6,  p.  163] 


J 


(4.31) 


is  a  quantity  of  some  interest.  It  represents  the  phase  difference 
between  the  observed  light  at  the  points  (u,,,  v2)  and  (u^,  v^). 
From  (4.32)  it  is  possible  to  show  that 

PA<Jr 


(A^)  = 


lf2  f(l  -  P2)1/2  +  e(lt  -  cos_1f?) 


(1  -  P2)3/2 


(4.35) 
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where 


R1  R2 

p  =  -=■  cos  A'Jr  +  -r  sin  A'Jr 

cr  ff 


(4.36) 


Finally,  the  mean  value  of  |Ai>'|  is  given  by  [Ref.  5,  p.  411  ] 


|  Air  |  =  n 


2  sin 


-1 


K  K 


-  u 


V 


2 

a 


(4.37) 
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V.  SPATIAL  AVERAGES 


The  preceding  chapter  considered  only  averages  over  an  ensemble  of 
scattering  surfaces.  However,  the  experimental  tools  that  an  observer 
must  use  are  spatial  averages  (i.e.,  spatial  integrations)  over  the 
observation  region  of  the  U-V  plane.  In  this  chapter  the  statistical 
properties  of  spatial  averages  are  examined,  with  particular  attention 
paid  to  discovering  conditions  under  which  spatial  averages  are,  with 
high  probability,  very  nearly  the  same  as  ensemble  averages. 

The  integral  of  primary  concern  here  is  an  average  of  the  random 
process  l(u,  v)  over  the  region 

-U  <  u  <  U,  -V  S  v  S  V  (5.1) 

of  the  observation  plane.  The  average  is  defined  by 

*  =  4UV  /v  jT  I<U-  V>  du  dV  <5-2) 

First,  the  expected  value  (over  the  ensemble  of  scattering  surfaces) 
of  the  spatial  average  is  examined  (5.2).  The  orders  of  expectation  and 
integration  may  be  interchanged  to  yield 

I  =  I(u,  v)  =  J a| 2  a2  (5.3) 

Thus  the  spatial  integration  yields,  on  the  average,  the  ensemble 
expectation  of  l(u,  v). 

It  is  pertinent  at  this  point  to  inquire  as  to  how  large  the  depar¬ 
tures  of  5  from  its  average  value  might  be  on  any  single  trial.  The 
dependence  of  these  departures  on  U,  V  and  on  the  statistical  structure 
of  l(u,  v)  is  of  particular  interest.  The  answer  to  this  question 
will,  of  course,  indicate  the  magnitude  of  the  difference  between  spatial 
averages  of  two  intensity  patterns  which  result  from  statistically 
independent  scattering  surfaces  (composed  of  similar  materials). 
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Therefore  the  variance  of  5  over  the  ensemble  of  scattering  surfaces 
is  calculated. 

To  begin,  note  that 


dui  du2  dVl  dV2 


(5.4) 


which  may  be  rewritten  as 


vg)  dux  du2  dvx  dv2 


(5.5) 


At  this  point  it  is  clear  that  the  variance  of  the  spatial  average  will 
depend  only  on  the  averaging  region  and  the  autocorrelation  function  of 
l(u,  v). 

Continuing,  (4.1l)  may  be  used  to  write 


'  W4  ’  *  “^2  '  “ 2 •  Vl  '  V2)|2d"‘  '“‘2  ',"1 


from  which  it  is  seen  that  the  variance  of  $  is  given  by 


|aj4  fVrrVf  2 

Var  U)  =  J-vJ' V,R«0<Ul  ’  V  Vl  ‘  ^  dU]L  d“2  ^  ^ 


A  change  of  variables  [Ref.  6,  p.  68]  allows  the  expression  to  be 


reduced  to 


lal4  r2V  r2V  ,  e  \  /  \  0 

U]  =  ~f0  f0  V  •  m)  (J  -  5v)  d*  d^  <5-8) 
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The  expression  (5.13)  contains  considerable  information.  One  con¬ 
clusion,  which  is  not  surprising,  is  the  following: 

As  the  size  of  the  integration  area  (uv)  increases,  the  spatial 
average  converges  to  the  ensemble  average  of  the  intensity  process. 

A  slightly  more  subtle,  but  related  conclusion  is: 

As  the  size  of  the  scattering  spot  (as  measured  by  y  and  y^) 

increases,  the  spatial  average  again  converges  to  the  ensemble 
average. 
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VI.  THE  DEPENDENCE  OF  SPARKLE  PATTERNS  ON  LIGHT 
FREQUENCY  CONTENT 


In  the  previous  chapters  ideally  monochromatic  light  sources  were 
assumed.  In  actuality,  the  output  of  a  laser,  regardless  of  its  type, 
consists  typically  of  a  number  of  spectral  components,  each  of  finite 
width.  Since  sparkle  patterns  were  not  observed  until  the  advent  of 
narrowband  laser  sources,  it  is  of  considerable  interest  to  determine 
the  exact  role  of  spectral  content  in  this  phenomenon.  This  chapter  is 
concerned  with  (l)  the  cross  correlation  between  sparkle  patterns  pro¬ 
duced  by  different  monochromatic  frequency  components;  (2)  some  statis¬ 
tical  properties  of  sparkle  patterns  produced  by  a  laser  spectrum 
consisting  of  a  number  of  monochromatic  components;  and  (3)  some  suf¬ 
ficient  conditions  for  observing  sparkle  patterns  when  the  spectral 
components  have  nonzero  width. 

A.  THE  CROSS  CORRELATION  OF  PATTERNS  PRODUCED  BY  DIFFERENT  FREQUENCY 

COMPONENTS 

Suppose  that  two  light  sources  with  different  frequencies  f^  and 
f2  cps  are  available.  Each  source  is  assumed  to  le  ideally  monochromatic, 
and  will  produce  a  sparkle  pattern  when  its  light  falls  on  a  diffuse 
surface.  It  is  pertinent  to  inquire  as  to  the  correlation  between  the 
two  patterns  at  each  point  (u,  v)  of  the  observation  region,  as  a 
function  of  the  frequency  difference  Af  defined  by 

Af  =  f2  -  t1  (6.1) 

assuming  that  in  both  cases  the  source  illuminates  the  same  area  of  a 
fixed  target  and  that  the  intensity  distributions  of  the  incident  spots 
are  identical.  If  the  correlation  is  high  for  all  (u,  v),  then  the 
two  sparkle  patterns  must  be  very  nearly  identical;  if  it  is  low,  the 
two  patterns  bear  little  resemblence. 
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From  the  definition  (4.15),  the  cross  correlation  of  the  intensity 
pattern  when  the  frequency  is  f  and  the  intensity  pattern  when  the 
frequency  is  fg  is  given  by 

P  (6u,  5v)  =  I(u,  v;  f  )  I(u  -  6u,  v  -  8v;  f  )  (6.2) 

112  1  2 


Bv  an  argument  similar  to  that  used  to  derive  (4.1l)  it  can  be  shown 
that 

ft,  (5u,  Sv)  =  |  a| 4  [o^o-g  +  |ft  (6u,  Sv ) 1 2_|  (6.3) 

12  L  12  J 


where 

F.e  (Su,  Sv)  =  e0(u<  v;  fx)  ej(»  -  Su,  V  -  6v;  fg)  (6.4) 

2 

while  0^/2  is  the  variance  of  the  real  and  imaginary  parts  of 

eQ(u,  v;  f^,  and  J^/2  is  the  corresponding  variance  for  eQ(u,  v;  fg). 

If  Af  is  sufficiently  small  to  assure  that  the  same  total  intensity  is 

2  2 

scattered  from  the  surface  in  both  cases,  then  cr^  and  are  identical, 

and  to  determine  ft_  it  suffices  to  determine  ft 

12  12 

To  continue,  what  is  the  correlation  between  the  pattern  l(u,  v;  f^) 
and  the  untranslated  pattern  l(u,  v;  fg)  at  each  point  (u,  v)?  It  is 
assumed,  for  simplicity,  that  the  observation  plane  U-V  is  parallel  to 
the  plane  defined  by  the  mean  surface  (X-Y),  although  the  laser  beam 
is  allowed  to  be  incident  on  the  surface  at  an  angle  $  with  the  normal 
to  the  X-Y  plane,  as  shown  in  Fig.  2.  The  Y  and  V  axes  are 
directed  upward  in  this  figure. 
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FIG.  2.  GEOMETRICAL  RELATIONSHIP  BETWEEN 
THE  LASER  BEAM,  THE  SCATTERING  SURFACE, 
AFD  THE  OBSERVATION  PLANE. 


In  Appendix  B,  the  desired  cross  correlation  is  found  to  be  given 
by  the  following  expression: 


(0, 

e12 


(6.5) 


where 


ftg  (8u,  8v)  =  eQ(u,  v;  f2)  e*(u  -  8u,  v  -  5v  j  fg)  (6.6) 


and  <t>z( 4)  is  the  characteristic  function  of  the  random  variables  zk, 
as  given  by  the  Fourier  transform  of  the  probability  density  function 
of  the  z^, 

/OO 

exp  (j2H|zk)  PzUk)  dzk  (6.7) 
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It  is  proposed  to  use  the  result  (6.5)  to  find  the  conditions  under 
which  the  two  spatial  patterns  £q(u,  v;  f^)  and  eQ(u,  v;  f2)  will  be 
highly  correlated.  For  high  correlation,  both  factors  in  (6.5)  must  be 
near  their  maximum  values.  The  magnitude  of  the:  characteristic  function, 
|$2(0(>  is  considerably  less  than  1^3  maximum  value  when  £  is  greater 
than  the  reciprocal  of  the  standard  deviation  crr.  of  p  (z  ).  Thus  a 
necessary  condition  for  high  correlation  is* 


Af  <  — 
°Z 


1 

cos 


(6.8) 


Turning  to  the  first  factor  in  (6.5),  it  can  be  shown  from  (4.5)  that 
|R  (ou,  Sv)j  is  considerably  less  than  its  maximum  when 

6u  >  or  &v  >  (6.9) 

f2  \  f2  "y 


where  t,  is  the  X-dimension  width  of  the  scattering  spot  and  l  is 
X  y 

the  Y-dimension  width.  Thus  a  necessary  condition  for  high  correlation 
is 


u  < 


Dc 

Af  l 

x 


and 


Dc 

aTT 

y 


(6.10) 


for  all  (u,  v)  included  in  the  observation  region. 

The  conclusions  of  this  section  can  be  summarized  as  follows: 

The  correlation  between  patterns  produced  by  different  frequency 
components  is  maximum  at  (u,  v)  =  (0,  0)  and  extends  over  a 
finite  region  of  the  U-V  plane.  The  magnitude  of  the  maximum 
correlation  is  determined  by  the  relationship  (6.8)  between  the 
frequency  difference  Af,  the  surface  roughness  c t,,  and  the 
angle  of  incidence  $>  of  the  laser  beam. 


This  result  is  quite  similar  to  the  so-called  Rayleigh  Roughness 
Criterion  [Ref.  4,  p.  411). 
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B.  STATISTICAL  PROPERTIES  OF  THE  SPARKLE  PATTERN  PRODUCED  BY  A  SUM 

OF  MONOCHROMATIC  LIGHT  COMPONENTS 

In  this  section  the  first-order  statistics  of  the  total  intensity 
pattern  are  found  for  two  special  cases  of  interest,  and  ar.  expression 
for  the  autocorrelation  function  of  the  total  intensity  pattern  is  also 
found . 

It  is  -ssumed  that  the  scattering  process  is  a  linear  one,  from 
which  it  follows  that  when  the  laser  output  consists  of  a  sum  of  mono¬ 
chromatic  frequency  components,  the  total  electric-field  pattern  is 
simply  the  sum  of  the  field  patterns  produced  by  the  individual  components 
separately.  The  total  electric  field  may  be  written,  from  (3.2),  as 

Hf  (u2  +  v2)] 

n'  _  / .  .  \ 


(6.11) 

where  the  f ^  are  the  frequencies  of  the  N  monochromatic  components 
involved. 

The  total  intensity  pattern  is,  from  the  definition  (3,5),  given  by 


N 

eT(u,  v;  t)  =  a  exp  (j2JIfntj  exp 
~  n=l  '  ' 


.V  .V 
,2  ^ 


IT(u,  v)  =  jaj  2^  2^  eXp 
n=l  Bal 


- 1  ;i 
n  n  1 


n(t  -  tj(u2  *  v2)' 


Dc 


«0(“.  V  fm> 


(6.12) 


If  the  intensity  detector  is  followed  by  a  low-pass  filter  sufficiently 
narrow  to  eliminate  all  but  the  dc  term,  (6.12)  becomes 

N  N 

IT<U’  v)=2l-!2  !fo(u>  V;  fn)i2=5^l(u’  v;  (6-13) 

n=l  n=l 
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Now  each  l(u,  v;  f  )  is  an  exponentially  distributed  random 
variable  which  may  or  may  not  be  significantly  correlated  with  other 
members  of  the  sum.  For  the  general  case  of  arbitrary  frequency 
spacings  between  spectral  components,  the  statistics  of  I^(u,  v)  are 
extremely  difficult  to  compute.  However,  two  limiting  cases  of  interest 
can  be  solved; 

1.  The  case  of  spectral  components  with  sufficiently  wide  frequency 
separations  to  produce  uncorrelated*  patterns:  and 

?  The  case  of  spectral  components  with  sufficiently  small  frequency 
separations  to  produce  perfectly  correlated  patterns. 

Suppose  first  that  the  minimum  separation  between  frequency  com¬ 
ponents  is  so  large  that  all  terms  of  (6.13)  are  uncorrelated  for  every 
(u.  v)  in  the  observation  region.  A  sufficient  condition  for  this  to 
be  true  is,  irom  (6.8), 

min  (f  -  f  }  »  4-  fl  +  — (6.14) 
/  n  m  o„  V  cos  V 
n^m  Z  '  ’ 


Further  suppose  that  the  scattered  light  intensities  resulting  from 
different  frequency  components  are  not  precisely  identical,  a  condition 
which  is  always  met  in  practice  due  to  the  fact  that  the  laser  spectral 
components  will  always  have  at  least  slightly  different  intensities. 
Then,  as  shown  in  Appendix  C,  the  probability  density  function  of  the 
total  intensity  I^,(u,  v)  is  given  by  a  weighted  sum  of  exponential 
functions, 


r 


-<q>  ■  ( 


E 

n=l 


i  2  2 

a  a 
n 


exp  - 


,2  2 
a  ff 
-1  n 


i^o 


otherwise 


(6.15) 


where  the  are  real,  constant  weighting  factors  which  are  evaluated 

in  Appendix  C. 


^Because  the  random  variables  l(u,  v;  fR)  are  composed  of  the  sum  of 

the  squares  of  two  normal  processes,  zero  correlation  implies  statistical 
independence. 
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Next,  suppose  that  the  maximum  separation  between  frequency  components 
is  so  small  that  all  terms  of  (6.13)  have  unity  correlation  in  the  obser¬ 
vation  region.  Sufficient  conditions  for  this  to  be  true  are 


max  (l 

n 

n,m 


f  )  «  — 

m  C7„ 


1  + 


COS  0 


-1 


(6.16a) 


and 


u  « 


Dc 


max  {f  -  f  ) 
n  m 


v  « 


Dc _ 

max  {f  -  f  )  & 


(6.16b) 


for  all  (u,  v)  in  the  observation  region.  Since  all  the  terms  of 
(6.13)  are  perfectly  correlated,  the  statistics  of  I^(u,  v)  are 
easily  seen  to  be  exponential,  with  a  probability  density  function 
given  by 


I  0  otherwise 

Finally  the  autocorrelation  function  of  the  total  intensity  pattern, 
as  given  by 


N  N 


5-t(6u,  5v)  =  I^,(u,  v)  I^,(u  -  6u,  v  -  ov)  -1SL,  Rj  (Su,  6v)  (6.18) 


n=l  m=l 


is  investigated.  From  (6.3)  it  follows  that 
N  N 


Rt(6u,  6v)  = 


=£]C  >-i4  Tv3™ +  k  (8u- 6v)i^| 

1  nm  J 

n=l  m=l  - 


(6.19) 
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while  from  Appendix  B  it  can  be  shown  that  if  (6.16b)  applies,  then 


Re  (feu,  6v)  =Re  (feu,  8v)  -2 —m)  (6.20) 

nm  m  —  \  / 


Thus, 

N  N 

Rt(5u,  Sv)  =  l-i4 

n=l  m=l 


2  2 
<J  CT 
n  m 


+  |P„J  IS,  (5u-  6v)|  (6.21) 


2 


where 


P  =  <t>„ 
nm  Z 


f  -  f 
n  m 

c  cos  © 


(6.22) 


is  independent  of  (feu,  6v).  Thus: 

the  autocorrelation  function  of  the  total  intensity  pattern  is 
given  by  a  weighted  sum  of  the  autocorrelation  functions  of  the 
individual  intensity  patterns. 


C.  SPARKLE  PATTERNS  PRODUCED  BY  SPECTRAL  COMPONENTS  OF  FINITE  WIDTH 

Suppose  that  the  laser  output  spectrum  consists  of  a  single  spectral 
component  of  finite  bandwidth  W,  as  shown  in  Fig.  3.  The  exact  shape 
of  the  spectrum  is  not  of  concern  here,  only  the  fact  that  it  extends 
over  a  range  W  cps.  In  this  section  it  is  not  required  that  a  low- 
pass  filter  eliminate  all  the  difference  frequency  terms,  but  rather 
the  observation  is  performed  directly  on  the  output  of  the  intensity 
detector. 

The  laser  output  field  strength  is  considered  to  be  a  narrowband 
random  process  in  time,  which  is  represented  as  an  amplitude-and  phase- 
modulated  sinusoid, 


n(t)  =  A( t)  exp  (j[u>0t  +  0( t) ]) 


(6.23) 
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where 


n(t)  =  strength  of  the  laser  output 

A(t)  =  random  amplitude  modulation,  of  bandwidth  W  (approxi¬ 
mately) 

6(t)  =  random  phase  modulation,  of  bandwidth  W  (approximately) 

o>  =  center  frequency  of  the  laser  spectrum. 

0 


OUTPUT 

SPECTRUM 


FIG.  3.  LASER  POWER  SPECTRUM  WITH  A  SINGLE 
COMPONENT. 


The  observer  is  assumed  to  examine  the  scattered  intensity  pattern  for 
a  time  T  sec.  This  observation  interval  is  subdivided,  for  analysis 
purposes,  into  a  number  of  subintervals,  each  of  duration  somewhat  less 
than  1/W  sec.  During  any  one  of  these  subintervals  the  laser  output 
is  practically  identical  to  a  monochromatic  sinusoidal  oscillation  of 
constant  amplitude  and  constant  frequency  o>,  given  by 

u>  =  u>q  +  9  (6,24) 

where  8  is  the  time  derivative  of  the  phase.  Both  the  amplitude  and 
the  frequency  of  the  oscillation  will  change  as  time  progresses  through 
the  subintervals.  Two  distinctly  different  cases  can  now  be  distinguished: 
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Case  1 


W  « 


cos  0 


(6.25a) 


Case  2 


W  »  — 

°z 


cos  0 


-1 


(6.25b) 


In  case  1,  the  maximum  excursion  of  the  oscillation  frequency  u)  is 
sufficiently  small  to  assure  that  the  intensity  patterns  produced  in  all 
subintervals  are  very  highly  correlated.  Thus  the  sparkle-pattern  struc¬ 
ture  does  not  change  significantly  from  subinterval  to  subinterval,  and 
an  observer  can  distinguish  a  stationary  sparkly  pattern,  regardless  of 
how  long  his  observation  time  T  may  be. 

In  case  2,  the  excursions  of  the  oscillation  frequency  are  so  great 
that  uncorrelated  patterns  may  be  produced  in  adjacent  subintervals,  and 
the  sparkle  pattern  will  "wash  out"  unless  the  observation  time  T  is 
less  than  l/W.  But  note  that  no  matter  how  wide  W  may  be,  for  a 
sufficiently  short  observation  time  a  distinct  sparkle  pattern  will  be 
observed. 

In  summary,  sufficient  conditions  for  observing  sparkle  patterns 

are 


W  «  — 

°Z 


_ l 

cos 


(regardless  of  T) 


(6.26a) 


or 


T  «  ~  (regardless  of  crz)  (6.26b) 

The  final  situation  to  be  examined  is  that  of  a  laser  spectrum  con¬ 
sisting  of  a  series  of  narrowband  spectral  components,  as  shown  in  Fig. 

4.  The  total  extent  of  the  spectrum  is  denoted  B  cps,  while  the 
width  of  any  one  of  the  spectral  components  is  approximately  W  cps. 
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OUTPUT 

SPECTRUM 


FIG. 


LASER  POWER  SPECTRUM  WITH  SEVERAL  COMPONENTS. 


In  this  case  the  laser  output  may  be  written  as  a  sum  of  narrowband 
processes, 


n(t)  =^Am(t)  exP  fml  +  ^(O!) 


(6.27) 


where  A^Ct)  and  0m(t)  are  bandlimiteri  to  approximately  W  cps. 
Reasoning  identical  to  that  used  for  a  single  output  component  shows 
that  the  following  conditions  are  sufficient  for  sparkle  patterns  to 
be  distinguished: 

L.  B  «  —  (l  +  1  -) 

crz  \  cos  $  / 

or 

2.  T  «  i 

In  the  present  case  it  is  not  sufficient  that 

w  «  (l  +  -  ) 

crz  \  cos  ^  ) 

for  under  this  condition  it  is  still  possible  for  different  spectral 
components  to  produce  uncorrelated  patterns  which  will  then  add  in 
proportions  that  ch.u.ge  from  subinterval  to  subinterval. 


(6.28a) 

(6.28b) 

(6.29) 
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APPENDIX  A.  FAR-FIELD  FOURIER-TRANSFORM  RELATIONS 


This  appendix  considers  the  relationship  between  the  field  distri¬ 
bution  E(x,  y;  t)  across  the  scattering  spot  and  the  field  distribution 
s(u,  v;  t)  in  the  observation  region.  The  observed  field  at  the  point 
(u,  v)  may  be  written  as  a  sum  of  contributions  from  the  incremental 
areas  AA  =  AxAy  on  the  spot.  As  AA  is  allowed  to  approach  zero,  the 
s.<m  passes  to  an  integral  which  may  be  written* 

/•<*>/•  E(x,  y;  t)  .  » 

l(u,  v;  t)  =  f  J  -  exp  ^  -  j  dx  dy  (A.l) 

where  r  is  the  distance  from  (x,  y)  to  (u,  v).  The  geometry  is 
shown  in  Fig.  S. 


FIG.  5.  THE  GEOMETRY  LEADING  TO  EQ.  (A.l). 


This  equation  follows  directly  from  Huygen's  principle. 
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Under  the  assumptions 

5«1'  D<<1 

for  all 

(x,  y) 

in  scattering  spot 

(A. 2a) 

|«1-  5«1 

for  all 

(«.  v) 

in  observation  region 

(A. 2b) 

(A.l)  lay  be  rewritten 

as 

exp  (j  <0  t)  />“/•  ,  n  x 

£(u,  v;  t)  = - j— -  J  J  Eq(x,  y)  exp  — J  dx  dy  (a.3) 


where  EQ(x,  y)  is  defined  by 

E(x,  y;  t)  =  Eq(x,  y)  exp  (j  cuQt)  (A.4) 

The  distance  r  is  seen  to  be  given  by 

*(t^)2}1/2  <*■»> 

The  binomial  expansion  may  now  be  used  to  expand  the  exponent  in  the 
integrand, 
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APPENDIX  B.  CROSS  CORRELATION  OF  FIELD  PATTERNS 
PRODUCED  BY  DIFFERENT  FREQUENCY  COMPONENTS 


This  appendix  investigates  the  cross-correlation  function  of  two 
complex  field  patterns  produced  by  different  monochromatic  frequency 
components.  Following  the  notation  of  Chapter  VI,  the  correlation 
function  of  interest  is  given  by 


As  the  result  of  the  statistical  independence  of  the  phases  of  different 
scatterers,  Eqs.  (B.l)  and  (B.2)  car.  be  combined  to  yield 


(B.3) 
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which  may  be  rewritten 


J_  Jsjt  )  s„(f  )  _ 

(8u.  Sv)  =  £  * - 5ii - -  exp  UlSk(*n>  -  Sk(f„ 


TIT 


k=l 


ifi(V  yk>l2  exp  {j  u  +  6u)  *k +  (-V-5  v  +  8v)  yy>]} 

(B.4) 


But  now  note  that 


r  2H(f  -  f  ) 
n  m 

rt  -  r 

n  _ m 

X  S\ 

Lj  ■ 

kCOS  ♦  i 

_ 

[cos  $  /J 

(B.5) 

where  *  is  the  angle  of  incidence  of  the  laser  beam  (see  Fig.  2),  z 

til  ” 

is  the  Z  coordinate  of  the  k  radiator,  and  4>  (|)  is  the  charac- 

Z 

teristic  function  of  the  random  variables  z^.  Finally,  if  the  frequency 
difference  f^  -  f^  is  sufficiently  small  to  assure  that  the  scattering 
cross  sections  of  the  elementary  scatterers  are  independent  of  frequency, 
comparison  of  (B.4)  and  (4.3)  shows  that 


'€  \  f 
m  \  m 


If  -  t 

ft.  (6u,  6v)=»  '-a  m 


where 


f  -  f 

u  +  6u,  -n  - — S  v  +  5v 


(B.6) 


K  — 


(6u;,  sv)  =  2 5f  lEi(xk’  yk)|2  exp  "dT  (xk6u'  yk5v')| 

k=l  L  J 


(B.7) 


is  the  autocorrelation  function  of  the  pattern  corresponding  to  f 
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APPENDIX  C.  FIRST-ORDER  STATISTICS  OF  THE  INTENSITY 
PATTERN  PRODUCED  BY  A  SUM  OF  WIDELY  SPACED 
MONOCHROMATIC  FREQUENCY  COMPONENTS 


In  this  appendix,  the  incident  light  beam  is  considered  to  be  com¬ 
posed  of  a  sum  of  monochromatic  spectral  components,  spaced  so  widely 
apart  that  the  intensity  patterns  produced  by  different  components  are 
statistically  independent.  As  indicated  by  Eqs.  (3.1l)  and  (6.13),  the 
total  intensity  pattern  is  simply  the  sum  of  a  number  oi  statistically 
independent,  exponentially  distributed  random  variables. 

It  is  a  well-known  result  of  statistics  that  the  characteristic 
function  of  a  sum  of  a  number  of  statistically  independent  random  varia¬ 
bles  is  simply  the  product  of  the  characteristic  functions  of  the  com¬ 
ponent  random  variables.  From  (3.1l),  the  n^  component  intensity 
pattern  has  a  probability  density  function  given  by 

r  /  \ 


I  ,  ,2  2  exM‘ 

p(!n)  =  <  ^  °n  \ 


,  ,2  2 

lal  V 


I  >  0 


(c.1) 


0 


otherwise 


The  corresponding  characteristic  function  is  [Ref.  7,  p.  221 } 


»n(s) 


(C.2) 


It  follows  that  the  characteristic  function  of  the  total  intensity 
pattern  is 


r<t)  *  n 


Ji  (i  -  J  2%12v:) 


(C.3) 
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where  N  is  the  number  or  frequency  components  involved.  Now  if  none 


ll 

I 


of  the  are  identical,  the  characteristic  function  has  only  single¬ 

order  poles,  and  a  partial-fraction  expansion  yields  the  simple  result 


N 


n=l 


1 


B 

n 


j  2K  |  a  j  2 


2 


I 


(C.4) 


where 


B  =  Um 

j|  •*  ^2fl 

An  inverse  Fourier  transform  of  (C.4)  yields  the  probability  density 
function  of  the  total  intensity  pattern, 


v _ 

B 

n  j 

f.  .  \ 

,  ,2  2  e  <P 

,  ,  .2  2 / 

a  cr 

n=l 

I  _i  n 

\  i  — '  n ' 

P  (IT>  =  ( 


(C.6) 


0 


otherwise 
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